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ABSTRACT 

We present new z ~ 8 galaxy candidates from a search over ~95 arcmin 2 of WFC3/IR data, tripling 
the previous search area for bright z ~ 8 galaxies. Our analysis uses newly acquired WFC3/IR imaging 
data from the CANDELS Multi-Cycle Treasury program over the GOODS South field. These new data 
are combined with existing deep optical ACS imaging to search for relatively bright (Mjjv < —19-5 
mag) z ~ 8 galaxy candidates using the Lyman Break technique. These new candidates are used to 
determine the bright end of the UV luminosity function (LF) of star- forming galaxies at z ~ 7.2 — 8.7, 
i.e. a cosmic age of 600 ± 80 Myr. To minimize contamination from lower redshift galaxies, we make 
full use of all optical ACS data and impose strict non-detection criteria based on an optical xl p t nux 
measurement. In the whole search area we identify 16 candidate z ~ 8 galaxies, spanning a magnitude 
range i^i60,AB = 25.7 — 27.9 mag. The new data show that the UV LF is a factor ~ 1.7 lower at 
Muv < —19.5 mag than determined from the HUDF09 and ERS data alone. Combining this new 
sample with the previous candidates from the HUDF09 and ERS data allows us to perform the most 
accurate measurement of the z ~ 8 UV LF yet. Schechter function fits to the combined data result 
in a best-fit characteristic magnitude of M*(z = 8) = —20.04 ± 0.46 mag. The faint-end slope is very 
steep, though quite uncertain, with a = — 2.06±0.32. A combination of wide area data with additional 
ultra-deep imaging will be required to significantly reduce the uncertainties on these parameters in 
the future. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: luminosity function, mass 
function 



1. INTRODUCTION 

Thanks to the unprecedented efficiency of the 
WFC3/IR camera on-board the Hubble Space Telescope 
(HST), the last two years have seen a remarkable 
progress in the exploration of galaxies within the first 
Gyr after the Big Bang. Immediately after its instal- 
lation, WFC3/IR has been used to obtain the deep- 
est NIR images ever seen as part of the HUDF0 9 pro- 
gram (PI: Illingworth; e.g. IBouwens" etaHl2011bh . This 
pushed the observational frontier of galaxies from z ~ 6 
well into the r eionization epoch, with first constraints 
even at z ~ 10 (IBouwens et al.ll20lTal : lOesch et al1l2012t 
iZheng et al.ll2012f ). 

In combina tion with deep ancil lary optical data from 
the HUDF dBeckwith et "all l2006h and GOODS surveys 
(jGiavalisco e tall l2QQ4D ,~the WFC3/IR data from the 
HUDF09 and the ERS (jWindhorst et alJl201lh programs 
allowed for the identification of more than 130 z ~ 
7 — 8 galaxy candidates to date using the Lyman Break 
Galaxy (LBG) selection technique (e.g. Ste idel et al.l 
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1996; Giavalisco 2002). Since the Lya absorption of 
the neutral inter-galactic medium shifts to > 1/im at 
z > 7, such galaxies can only be detected in the NIR 
data and are completely invisible in the optical. For ex- 
ample, z ~ 8 galaxies can be selected as Yios-dropouts, 
based on their very red Y105 — J125 colors and optical 
non-detections. 

The first z ~ 7 — 8 galaxy candidates have 
already been used to measure th e UV LF at 
these r edshift s (e.g. [O esch et al. 2010c; Bouwens et al 



2010a, 201l fl iMcLure et all [201G : I Bunker et all 
Fin kelstein et all 120101 : lYan et alJ l2010t Wilkins 



2MQl : lLorenzoni eT^^oTTI : iBradlev et all I2MI ). Ad- 
ditionally, these candidates allowed for first estimates 
of the physi cal parameters of z > 6 galaxies, such as 
their sizes (|Oesch et al. 2010 a]), UV continuum slopes 
(e.g. Bouwe mTeTaTI l2010bL l2012bl : iFinkelstein et al.l 
I2010L IgQllt iDunlop et all |2012[ ) . rest-frame optical col- 



et all 
20T0l: 



et al.l 



ors (Gonzalez et al. 2011^ a nd eve n ages and stel- 

lar masses (e.g. Labbe e t all l2010bllal; [Gonzalez et al ' 



2010, 2011; Finkelstein et al. 2010; Schaerer & d e Barrosl 
2010). Furthermore, z ~ 8 proto-cluster candidates 



have been id entified in pure parallel WFC3/IR data 
(jTrenti et al.ll2012h . 

While the first deep WFC3/1R data sets provide good 
constraints on the faint population at z > 7, they only 
probe a limited volume and therefore result in poor 
constraints on the much less abundant bright galaxies 
around the exponential cut-off of the LF. At z ~ 7, 
constraints on this bright cut-off can be obtained us- 
ing ground-b ased, wide-area data from, e.g., Subaru or 
the VLT (e.g. IQuchi et al.ll2009l : ICastellano et al.ll2010d : 
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Mickey et al.ll2Q10h . However, beyond z ~ 7, the galaxy 
population is too faint, and the sky is too bright to effi- 
ciently detect z > 8 galaxies from the ground with cur- 
rent facili ties, resulting only in u pper limits at Mjjv < 
-22 mag (|Castellano et al.ll2Q10bh in the field. Note that 
ground-based searches behind lensing clust ers resulted in 
a few potential z > 8 candidates (see e.g. lLaporte et al.l 
1201 lh . and in limits on the bright end of the z ~ 9 LF 
(jLaporte et al.ll2Q12f ). 

Due to the limited volume probed at z ~ 8 with 
current deep HST data, the bright end of the z ~ 8 
LF has thus remained quite uncertain, subject to large 
cosm ic variance. Fo r exam ple, in the HUDF09-2 field 
data Bo uwens et al.l (j2011b[ ) identified t wo very bright 
(iJifio ab 21 26 mag) sou rces (see also Wil kins et al.l 
I2011at iMcLure et al.l l2011h . which possibly biased the 
previous measurement of the z ~ 8 UV LF towards some- 
what higher values. It is therefore essential to study 
larger area WFC3/IR data to quantify how represen- 
tative these bright candidates are of the z ~ 8 galaxy 
population. 

There are a few ongoing and planned HST programs 
that can be used for this task. First, the Multi-Cycle 
Treasury program CANDEL S is obtaining "^105 data 
over the two GOODS fields (jGiavalisco et al.ll2QQ4D . re- 
sulting in a total search area for z ~ 8 galaxies of 
~ 300 arcmin 2 . Second , the pure parallel programs 
BOR G (Flrenti et al.llgoTTHlTrenti et alll2012h and HIP- 
PIES (jYan et al.l l20 111 ), are designed to find rare bright 
z ~ 8 galaxy candidates, sampling a combined total of 
about 400 arcmin 2 to varying depths. Before the ad- 
vent of JWST, the combination of different WFC3/IR 
data sets probing different scales will likely be the only 
way to obtain a we ll-sampled UV LF at z ~ 8 (see e.g. 
iBradley et aTll2012D . 

In this paper, we take advantage of the full Yios-band 
data obtained over the GOODS-South field as part of the 
CANDELS program to identify relatively bright z ~ 8 
galaxy candidates at absolute magnitudes Mjjv < —19-5 
mag. We combine these with our previous z ~ 8 searches 
in the ultra-deep data from the HUDF09 as well as the 
ERS programs to significantly improve the sampling of 
the bright end of the z ~ 8 UV LF. With the inclusion 
of the CANDELS field, we essentially triple the search 
area for bright z ~ 8 galaxies in GOODS-South relative 
to our previous analysis, resulting in much more reliable 
constraints on the shape of the bright end of the LF 
and the luminosity at which an exponential cutoff might 
occur. 



An accurate measurement of the bright part of 
the UV LF is critical for enabling meaningful com- 
parisons with galaxy evolut i on mo d els at these 



redshifts (e.g. Finlator et al.l l2011bt iFinlator et al. 


2011; 


Trenti et al.l 120101: lLacev et al.l 120111: Jaacks et al. 


2012; 


Daval & Ferrara 2012: For ero- Romero et al.l 120101: 


Salva 


erraetal. 2011; Muhoz 2012). The bright end 



potentially contains important information on the star- 
formation efficiencies, formation timescales and duty cy- 
cles of galaxies within dark matter halos in the first few 
hundred Myr of cosmic time. 

Furthermore, a well sampled bright end of the LF helps 
to break degene racies in fitting the Schechter function 
( Schechterl ll976[ ) . Even if shallower data sets do not pro- 
vide real information on the faint-end slopes, they can 



still decrease the uncertainties on the faint-end slope, 
assuming that the shape of the LF can accurately be 
de scribed by a Schech ter function (e.g., see discussion 
in lBouwens et al.1l2008f ). Accurate estimates of the faint- 
end slope are critical for obtaining realistic constraints on 
the number of ionizing photons emitted by the ultra- faint 
galaxy population at these redshifts. Within the current 
uncertainties it is unclear whether the faint galaxy popu- 
lation was lu minous enough to reionize the universe alone 
at a > 7 (e.g. Bolton fc Haehne lt 2007; Oes ch et al.l20Q9t 
Robertson et al.ll2010l : iBouwens et al.ll2Q12at iShull et all 
20lTh . 

Finally, the identification of brighter z ~ 8 galaxy can- 
didates is crucially important to provide candidates for 
spectroscopic follow-up. Spectroscopy has proven to be 
extremely challenging, mainly due to the intrinsic faint- 
ness of the sources, the abundance of bright night-sky 
lines, and the absorption of a significant fraction of the 
Lya flux of these galaxies by the neutral inter-galactic 
medium (IGM). Therefore the progress in spectroscopic 
follow-up has been relatively slow, with only a handful 
of confirmed sources at z ~ 7 (see e.g. ISchenker et al.l 
12012b iPentericci et al.ll2QlH lOno et al.lT2012D . At z ~ 8, 
possibly the best chance for spectroscopic confirmation of 
the redshifts of candidates is with upcoming multi-object 
NIR spectrographs targeting simultaneously a few bright 
sources identified in contiguous wider area imaging. 

This paper is organized as follows. After the presenta- 
tion of the data in Section [2] we explain our source selec- 
tion criteria and present the z ~ 8 candidates in Section 
[3l These are then used to derive constraints on the bright 
end of the z ~ 8 UV LF in Section HI We will refer to the 
HST filters F435W, F606W, F775W, F814W, F850LP, 
F105W, F125W, F160W as £435, ^606, *775, 
^105? J125? #160, respectively. Throughout this paper, we 
adopt Vt M 0.3, Q A = 0.7, Hp = 70 km s ^Mpc" 1 , i.e. 
h = 0.7 (WMAP-7: lKomatsu^"aIll2011l| . Magnitudes 
are given in the AB system (jOke fc Gunrifl983l ). 

2. THE DATA 
2.1. WFC3/IR Data over GOODS-South 

In this paper, we analyze the complete public 
WFC3/IR data obtained over the GOODS South 
field as part of the Mult i- Cycle Treasury program 
CANDELS (PI: Fabe r /Ferguson; iGrogin et all I20TTI : 
iKoekemoer et all l2011f ) and then combine it with the 
deeper HUDF09 and ERS datasets. The completion of 
the Y105 imaging over the CANDELS GOODS South 
field (see Figure [1]) has made this a unique data set for 
studying the bright end of the z ~ 8 UV LF based on a 
Yios-dropout selection. 

The CANDELS data is split in two parts, CANDELS- 
Deep and CANDELS- Wide. The first covers the cen- 
tral part of GOODS South in 3 x 5 tiles with ~3.5 
orbits in each filter Y105, J125 and Hiqq. These data 
cover ~ 65 arcmin 2 , reaching down to Hiqo^ab = 27.7 
mag. The CANDELS- Wide program consists of data in 
the same three WFC3/IR filters as for the CANDELS- 
Deep. These were fully acquired already by end of 
March, 2011 and cover 8 WFC3/IR pointings (~ 35 
arcmin 2 ) to ~ 1 orbit depth in each filter, thus reach- 
ing to H 160iAB = 27.0 mag. 

In addition, we combined all the public WFC3/IR data 
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TABLE 1 

5cr Depth and Area of Data Used in this Paper 



Field Area [arcmin 2 ]"^ B 435 V606 ^775 ^814* ^850 ^105 ^125 #160 

CANDELS-Deep 60.2 27.9 28.1 27.6 28.5 27.6 27.9 28.0 27.7 

CANDELS-Wide 34.2 27.9 28.1 27.6 28.0 27.6 27.2 27.3 27.0 



Note. — The limits correspond to 5<r variations in the sky flux measured in a circular aperture 
of 0"25 radius, i.e. no correction to total magnitude was performed. The full CANDELS-Deep 
is about 3.5 orbits per filter, and CANDELS-Wide one orbit per filter, leading to a difference of 
about 0.7 mag. 

t The quoted area corresponds to the part of the image where imaging in all three WFC3/IR 
filters is available as well as ACS in S435, V606, 2775, £850 . 

* ^814 is used for confirming the optical non-detection of candidates using a 2a non-detection 
criterion. All candidates listed in table [2] do have /si4 coverage. 




Fig. 1. The WFC3/IR F105W data over the GOODS South 
field that are used in this analysis. These data include all 
WFC3/IR data taken over this area before the end of 2011. The 
different fields are color coded by exposure time, with deeper ar- 
eas being darker. The CANDELS-Deep covers ~ 65 arcmin 2 with 
~3.5 orbits of F105W data. The outline indicates the area of the 
optical ACS coverage of GOODS. The CANDELS-Wide field spans 
34 arcmin 2 and is only of 1 orbit depth. The central part covered 
by the HUDF i s omitted from our prim ary analysis since it was an- 
alyzed fully by Bouwens et al. (2011b). However, we do combine 
the CANDELS z ~ 8 candidates with the ones w e previously se- 
lected in the deep HUDF09 and the ERS fields (see Bouwens et al. 
2011b) when constraining the z ~ 8 LF. For more information on 
the data used here see Section [2] and Table [T] 

that have been taken over these fields before March 2012 
from other programs. In particular, we included the 
imaging data of the supernovae follow-up program of 
CANDELS (PI: Riess), which adds imaging over a few 
pointings of CANDELS Deep. The final exposure map 
of all the Y105 data included is shown in Figure [TJ The 
central part covered by the ultra-deep HUDF09 program 
was omitted from our analysis of the CANDELS data 
in order not to replicate t he bright candidates o f that 
field given in the study of iBouwens et al.l (j2011bf ). We 
will use those candidates , how ever, when com puti ng the 
total z ~ 8 LF in section l4~4l (see also section [22} • 

The WFC3/IR data was reduced using standard tech- 
niques. The pipeline-processed science frames were ob- 
tained from MAST, and were subsequently registered to 



the existing optical ACS data at a pixel scale of 0"06. All 
input images were inspected visually for satellite trails 
and other artefacts, such as loss of the guide star. Pix- 
els affected by persistence from previous observations, 
as identified from the persistence masks provided by 
STSdQ, have been flagged and removed from our reduc- 
tions. For the last few visits of F105W data, we found 
that the cosmic ray subtraction of the archived data was 
not satisfactory. We therefore masked all pixels that were 
affected by cosmic rays in the archived data which re- 
sulted in much cleaner reductions. 

The optical ACS data used here include the original 
GOODS optical imaging as well as additional data avail- 
able over these fields from supernova follow-up programs. 
These reach about 0.1 — 0.3 mag deeper in zs^o than the 
v2.0 reductions of GOODS. This is critical for excluding 
interlopers to the z ~ 8 galaxy selections. We restricted 
our analysis to the part of the WFC3/IR data where 
imaging in all three filters (as well as full ACS cover- 
age) was available. This results in 60.2 arcmin 2 and 34.2 
arcmin 2 for CANDELS Deep and Wide, respectively. 

Finally, we also include a full reduction of all the ACS 
F814W filter data available over GOODS South. Such 
data was mainly taken in parallel to other observations of 
several programs. In particular, most of the CANDELS 
WFC3/IR observations obtained F814W data in parallel. 
However, over the years, several such programs have been 
conducted, and the resulting data reach to a non-uniform 
depth of igi4 ~ 28.0 — 28.5 mag, i.e. they are deeper than 
the original GOODS 2775 data, and in several parts of 
the field are essentially the deepest filter data now. The 
specifics of all the data used in this paper are summarized 
in Table [Q 

2.2. Deeper WFC3/IR Data 

For accurately constraining the UV LF at z ~ 8 it 
is crucial to probe the galaxy population over as large a 
dynamic range in luminosity as possible. In the later sec- 
tions of this paper, when we compute the LF, we there- 
fore include all the z ~ 8 galaxy candidate s that we previ- 
ously identified in Bouwens et all (|2011bf ). These candi- 
dates were selected from the three ultra-deep WFC3/IR 
pointings of the HUDF09 program (PI: Illingworth) as 
well as from ~40 arc min 2 of WFC3/IR Ear ly Release 
Science data (ERS; see Wind horst et al.ll2011f ). In total, 
these are 59 z ~ 8 galaxy candidates with magnitudes in 

8 http:/ /archive. stsci.edu/prepds/persist/ 
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Fig. 2. — L eft - The color sele ction criterion used to identify z ~ 8 galaxies. The color selection (blue shaded area) is identical to the 
one used by Bouwens et al. (2011b). This ensures that the resulting galaxy samples can be readily combined in the analysis. Additionally, 
we impose strict non-detection requirements based on the optical x 2 values (see text). This enables us to reliably select galaxies to quite 
blue Y105 — J125 colors. The redshifted color tracks of different types of galaxies are shown as solid (star-forming) and dashed lines (earlier 
types). The latter are shown up to z = 3, while the tracks of star-forming galaxies are extended to z > 8. Their locations at z = 7, 7.5, and 
8 are labelled and marked with small dots. The three blue lines correspond to different amounts of dust reddening E(B — V) = 0, 0.1, 0.2 
(left to right). The expected location of ultra-cool M, L, and T dwarf stars is indicated by a green hatched region. The colors of the z ~ 8 
candidates identified in this paper are shown as black squares with the usual la error bars. Right - The optical % 2 for galaxies with 
^160, AB — 27.0 mag against J125 — Hiqq color. This constraint is particularly effective in excluding intermediate redshift galaxies (see 
Section [331 and Figure [3]). The lines correspond to the same galaxy types and redshift ranges as in the left panel. The measurements for the 
galaxy candidates identified in the CANDELS data are shown as black squares. Upper limits are shown for galaxies with log Xopt < —1.2. 
The size of the plot symbols represent the measured Y105 — J125 colors. Only one source with Y105 — J125 < 0.9 lies near to the track of 
passive, intermediate redshift galaxies in this Xopt vs - color diagram. 



the range i^i60,AB = 26.0—29.4 mag over the total area of 
~ 53 arcmin 2 . For more information on t hese candidates 
and their selection, we refer the reader to Bouwens et al. 
(|2011bD . As we describe belo w, we use essentially iden- 
tical selection procedures to iBouwens et al.l (j2011bf ). to 
ensure that we can combine the new LF results and the 
older ones without corrections or biases. 

3. SOURCE SELECTION 
3.1. Catalog Construction 

Sour ce catalogs are obtained with the SExtractor pro- 
gram (jBertin fc Arnoutsl 119961 ). which is used to de- 
tect galaxies in t he square-root of chi-square image 
(jSzalay et al.l Il999h computed from the J125 and Hiqq 
data. We perform matched aperture photometry on 
point spread function (PSF)-matched images. The col- 
ors used here are based on small elliptical apertures (1.2 
Kron; these are typically nearly round and about 0.2" 
long axis), and total magnitudes are measured in stan- 
dard 2.5 Kron apertures (typically 0.4" radius), corrected 
to total fluxes using the PSF encircled energy of an equiv- 
alent aperture to account for flux loss in PSF wings (typ- 
ically ~ 0.15 — 0.2 mag). 

The input RMS maps were scaled to properly represent 
the flux variations in 0'/25 radius apertures, which were 
used to establish the detection significance of sources. 
The RMS scaling was done based on the detected vari- 
ation in 1000 random apertures per WFC3/IR frame on 
empty sky regions after 3<r clipping. This procedure 
ensures that the weight maps correctly reproduce the 
actual noise in the images. Subsequently only sources 
with signal-to- noise ratios (S/N) larger than 4.5 in Hiqq 



are considered (as measured in 0'/25 radius circular aper- 
tures). 

3.2. z ~ 8 Color-Color Selection 

Thanks to the strong IGM absorption in the rest-frame 
UV, z ~ 8 galaxies can be selected in broad-band color- 
color diagrams. At z > 7 the Lya absorption shifts into 
the Y105 band, rendering galaxies red in Y105 — J125, while 
the J125 — #160 color of star-forming galaxies is still ex- 
pected to be blue. The selection is discussed here, but 
the most challenging aspect is to minimize the contami- 
nation from low redshift objects, as we discuss in Sections 
[Q and 1331 

We adopt th e same color selection cr iteria for our z ~ 8 
selection as in IBouwens et al.l (j2011bl ) , which will allow 
us to directly combine the previous data with the can- 
didates identified here. Specifically, the selection criteria 
are: 

(I105- ^125) > 0.45 
(J125 -#i6o)<0.5 

Additionally, we require sources to be detected at 4.5a in 
both J125 and iifi6o> and at 5<r in at least one of these two 
bands. Finally, we require candidates to be detected at 
< 2<T in all the optical ACS data, and also require candi- 
dates to meet an optical xl p t A ux constraint. Both tests 
play a major role in removing potential low-redshift con- 
taminants (see Section l3^3|) . Finally, we checked the isi4 
images of all sources and additionally required galaxies 
to be 2(7 non-detections also in Isia- This last test helped 
to eliminate three sources which are likely contaminants. 
The above color criterion is illustrated in Figure [2j 
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where we show the expected colors of star-forming galax- 
ies using 100 Myr old st ellar population models from 
iBruzual fc Charlotl (|2003l ). These are additionally red- 
dened with small am ounts of dust using the dust law of 
iCalzetti et all (12000). The IGM absorption is modeled 
using iMadaul (j!995f ). Our color criteria start to select 
such galaxies z > 7.2. The selection then peaks around 
z ~ 8 and falls off due to the combination of an increased 
distance modulus and flux reduction as the Lyman Break 
approaches t he filter bandpass lim it at higher redshift 
(see Fig. 4 in Bouwe ns et al.|[2011bh . The mean redshift 
of our sample is (z) = 7.9, with 80% of galaxies expected 
to lie at z = 7.2 - 8.7 (see Fig [5]). 

In Figure [2l we also show t he co lors of lower redshift 
galaxies using Colem an" et al.l (119801) tem plates, as well as 
cool dwarf stars (|Burgasser et al.l l2004), which could po- 
tentially contaminate our selection. Despite some over- 
lap in this color-color diagram with z > 7 galaxy can- 
didates, we expect the contamination of such sources to 
be small due to our strict optical non-dete ctio n criteria 
(see next section and discussion in section [33]) . This is 
illustrated in the right panel of Figure [2j where we show 
the expected x%t values (defined in the next section) as 
a function of J125 — Hiqq color for different galaxy types. 
Most intermediate redshift galaxies are expected to be 
well detected in the optical data, resulting in log xl p t > 1- 
Only completely quiescent galaxies would lie close to our 
selection box. The figure also includes the measure d val - 
ues of the CANDELS z ~ 8 candidates (see section [33]) . 
The sizes of the symbols represent their Y105 — J125 colors. 
Our selection includes only three galaxies with relatively 
blue Y105 — J125 colors with positive x 2 op t values between 
1 and 3. These would have the highest probability of 
being lower redshift contaminants. 

As a further step against interlopers, we check 
each candidate for bright detections in the Spitzer 
IRAC 3. 6 ^m images available over GOODS-South (e.g. 
iDickinson et al.l [2003) . This has proven to be a very 
effective test for removing dusty contaminants in very 
high-z galaxy searches in which only limited information 
on the UV- continuum slope of galaxies is available (see 
e.g. lOesch et al.l [2012). Young, star- forming galaxies at 
z > 7 with small amounts of dust reddenin g are expected 
to sh ow colors Hiqo — [3.6] < 1 (see also iGonzalez et aTl 
|2011[ ). Any galaxy with substantially redder colors is 
therefore likely to be a lower redshift contaminant. 

Indeed, in constructing our sample, we identified one 
bright (#160 ab — 24.3 mag) edge-on spiral galaxy (at 
03:32:14.72,-27:46:21.6), which satisfied our WFC3/IR 
color and optical non-detection criteria. However, this 
source is extremely luminous in the IRAC data, with 
#160 — [3.6] = 2.3, and is also clearly detected even in 
the [8.0] band, as expected for a heavily dust obscured 
source at intermediate redshift. We therefore removed it 
from the potential z ~ 8 galaxy sample. All other sources 
where we could measure IRAC fluxes (i.e. which were not 
completely blended with bright foreground sources) were 
consistent with the limit #160 — [3.6] < 1. 

3.3. Optical Non- detections Using x 2 op t 

As indicated above, one of the main challenges in se- 
lecting robust LBGs is to remove intermediate redshift 
contaminants. Given the high efficiency of WFC3/IR the 
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Fig. 3. — The optical x 2 non-detection criterion in the CAN- 
DELS Deep data. The histograms show the different distribution 
functions for empty sky positions (gray filled), as well as what is 
expected for contaminants based on our simulations. True z ~ 8 
sources will have a distribution like the gray empty sky distribu- 
tion. The orange distribution is based on simulations using galaxies 
from the HUDF09 to which we applied Gaussian flux scatter ap- 
propriate to our ~ 2 mag shallower CANDELS data. The dark red 
histogram corresponds to simulations using brighter galaxies in the 
CANDELS field itself, which are dimmed and have the appropri- 
ate flux scatter applied. The distributions from the two simulations 
are very similar. The dark red filled histogram represents the frac- 
tion of all contaminants that satisfy the color and non-detection 
criteria, but which lie below the adopted Xopt nrQ it of 3.0 (see 
text), and which would be confused with real z ~ 8 objects. With 
this cut we are able to reduce the contamination rate by a factor 
~ 3-4 (the ratio of the dark red area to the total). The residual 
contamination in our fields is expected to be about 1 source in the 
CANDELS Deep data and < 0.1 in CANDELS Wide giving a total 
contamination rate of about 10%. 

available ancillary optical data in most of the WFC3/IR 
search fields is not appreciably deeper than the new IR 
data (see e.g. Table [1]). Thus at the faintest magnitudes 
even a 2<j-nondetection criterion, which is the standard 
that is used in LBG selections, is not sufficiently effective 

to eliminate faint, dusty interlopers. 

In previous papers (e.g. lBouwens et a"T1l2011bl ), we have 
developed an efficient method for eliminating \ow-z con- 
taminants by making full use of all the information in 
the optical data. In particular, for each galaxy we com- 
pute an optical pseudo xl P t value from its aperture flux 
measurements of all optical bands as 



Xopt 



£SGN(/i)(/iM 



where i runs over .8435 1^06 2775 and ^850, and SGN 
is the sign function, i.e. SGN(x) = — 1 if x < and 
SGN(x) = 1 if x > 0. 

For real high-redshift candidates the Xopt distribution 
is expected to be centered around zero, while for contam- 
inants the distribution is skewed toward positive values. 
Therefore, by removing galaxies with values above some 
limiting value Xum ^ * s possible to significantly reduce 
the number of contaminants in LBG samples. 

The adopted limiting value is derived from two differ- 
ent sets of simulations. The results are shown in Figure 
[3] where we compare the observed x 2 op t t° those from our 
simulations to establish the limiting value. 

The first simulation is based on the HUDF09 data, 
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TABLE 2 

Table of z ~ 8 Galaxy Candidates in the CANDELS GOODS-South Data* 



ID a S H i 6 o J125 - #160 Vios - ^125 S/N (if 160 7^125) Alternate ID a 



CANDELS-Deep 



CANDY-2499448181 


03:32:49.94 


-27:48:18.1 


25.70 


± 


0.09 


-0.0 ±0.1 


1.0 


± 


0.1 


22.0/30.0 


ISCL085 


CANDY-2402644099 


03:32:40.26 


-27:44:09.9 


26.19 


± 


0.16 


-0.2 ±0.3 


0.5 


± 


0.3 


4.6/6.4 




CANDY-2468950074 


03:32:46.89 


-27:50:07.4 


26.20 


± 


0.14 


-0.0 ±0.1 


0.5 


± 


0.2 


11.0/12.2 




CANDY-2209751370 


03:32:20.97 


-27:51:37.0 


26.46 


± 


0.14 


-0.2 ±0.2 


1.2 


± 


0.3 


7.8/8.9 


AUTO_212 


CANDY-2272447364 


03:32:27.24 


-27:47:36.4 


26.70 


± 


0.20 


-0.2 ±0.2 


1.0 


± 


0.3 


6.1/8.3 




CANDY-2320345371 


03:32:32.03 


-27:45:37.1 


26.74 


± 


0.26 


-0.0 ±0.2 


0.7 


± 


0.2 


13.1/13.0 




CANDY-2350049216 


03:32:35.00 


-27:49:21.6 


26.90 


± 


0.16 


0.1 ±0.2 


> 


2 


.3 


7.1/7.7 


ISO_063 


CANDY-2139147577 


03:32:13.91 


-27:47:57.7 


27.09 


± 


0.28 


0.1 ±0.3 


0.6 


± 


0.4 


6.2/5.0 




CANDY-2243349150 


03:32:24.33 


-27:49:15.0 


27.11 


± 


0.20 


-0.3 ±0.2 


0.5 


± 


0.4 


5.9/9.2 




CANDY-2181952456 


03:32:18.19 


-27:52:45.6 


27.14 


± 


0.19 


-0.2 ±0.2 


1.1 


± 


0.4 


7.0/8.9 


AUTO_204 


CANDY-2209848535 


03:32:20.98 


-27:48:53.5 


27.22 


± 


0.28 


0.1 ±0.2 


1.3 


± 


0.5 


8.2/8.4 


ISCL071 


CANDY-2209246371 


03:32:20.92 


-27:46:37.1 


27.31 


± 


0.25 


0.1 ±0.3 


0.8 


± 


0.6 


5.8/6.5 




CANDY-2277945141 


03:32:27.79 


-27:45:14.1 


27.76 


± 


0.35 


-0.3 ±0.3 


1.6 


± 


0.6 


4.5/6.5 




CANDY-2432246169 


03:32:43.22 


-27:46:16.9 


27.88 


± 


0.43 


0.2 ±0.5 


0.9 


± 


1.0 


4.5/5.0 





CANDELS-Wide 

CANDY-2379552208 03:32:37.95 -27:52:20.8 26.45 ±0.12 0.5 ± 0.2 1.1 ± 1.0 9.7/7.7 

CANDY-2408551569 03:32:40.85 -27:51:56.9 27.30 ± 0.22 -0.3 ± 0.3 1.4 ±1.0 5.4/8.8 

* Limits are la. 

a IDs refer to Yan et al., ApJ submitted, eprint arXiv:1112.6406v2 (see also appendix). 



where the available optical and WFC3/IR data are ~ 
1.5 — 2 mag deeper than in the fields studied here. At a 
given Hiqo^ab magnitude, this allows us to check what 
fraction of sources would contaminate our sample, if they 
were observed at the shallower depth of our data. We use 
a Monte-Carlo simulation in which we apply Gaussian 
scatter to the fluxes of the HUDF09 sources scaled to 
the depth of the CANDELS Deep and Wide data. We 
then compute the xl p t values for all sources which did not 
satisfy our selection criteria in the ultra-deep data, but 
which would have been selected and would have passed 
the optical 2a non-detection criterion. The distribution 
of Xopt values of these contaminants is shown in Figure 
[3] as an orange histogram. 

For the second simulation, we use brighter sources di- 
rectly from the CANDELS fields which are dimmed to 
fainter magnitudes. Thanks to the wider area relative to 
the HUDF, this allows us to probe rarer interloper pop- 
ulations and so is complementary to the HUDF09 sim- 
ulation. We only use sources with i^i60,A# = 24 — 25.5 
mag as input, in order to keep the population as close 
as possible to the magnitude range of our candidates. 
Again, for each source satisfying the color selection after 
applying Gaussian noise, we compute the value of Xopt- 
This is shown in Figure [3] as a dark red line. 

As is apparent from the Figure, the Xopt distribution 
of interlopers is significantly skewed to positive values 
relative to regions of empty sky as expected for real z ~ 8 
sources. After some consideration we adopted a limit 
Xopt < 3.0 as being a good discriminator between z ~ 8 
galaxies and contaminating sources. By using this limit 
of Xopt < 3.0, it is possible to reduce the contamination 
rate by a factor of ~ 3 — 4 relative to a simple 2<j non- 
detec tion criterion in all opt ical bands (see also Appendix 
D of Bouwens et al.ll2011bf ). At the same time, we lose 
just ~15% of real z ~ 8 sources. We correct for this 
small loss in the estimate of the selection volume. 

The expected residual number of contaminants, after 



applying the xl p t limit, is 1.0 ± 0.5 sources in the CAN- 
DELS Deep data, and < 0.1 in the whole CANDELS 
Wide field. Based on these simulations, we thus esti- 
mate a contamination rate that is only ~ 10% for our 
full sample. Note that, given our ignorance of the LF 
of interlopers and cosmic variance in the input galaxy 
sample for our photometric scatter simulations, it is con- 
ceivable that the true contamination rate is as high as 
15 - 20%. 

3.4. The Bright CANDELS z~8 Candidates 

Applying all the selection criteria outlined above, we 
identify 16 new z ~ 8 galaxy candidates in the full 
CANDELS Y105 data over the GOODS-South field; 14 
in CANDELS Deep, and two in CANDELS Wide. Their 
properties are listed in Table [2j and individual images for 
all candidates in the WFC3/IR bands as well as a stack 
of all the optical ACS images are shown in Figure HI 

The candidates span a range of i?i60,AB = 25.7 — 27.9 
mag. Our brightest source is thus among the most lumi- 
nous z ~ 8 candidates known to date. It is ~ 0.5 mag 
brighter than all other CANDELS sources, and would be 
an ideal target for future spectroscopic follow-up. 

Such a luminous object raises an interesting possibil- 
ity. If we assume that the UV luminosity of galaxies 
scales with halo mass, such a bright z ~ 8 galaxy is 
expected to lie in an over-dense environment with ~ 5 
fainter galaxies w i thin a diameter of 60 arcsec (see e.g. 
iTrenti et a,l-1l2012t iMuhoz fc Loebll2008D . Unfortunately, 
the source lies close to an extremely bright star, reduc- 
ing the detectability of fainter companions significantly. 
Nevertheless, we would have expected to see a few fainter 
companions around this source. None are found, how- 
ever. Unfortunately, the proximity to the bright star also 
makes it impossible to obtain IRAC flux measurements 
from the current Spitzer data to strengthen or refute the 
high redshift solution of this source. Deeper optical or IR 
data (or a spectrum) may thus be needed to fully resolve 
the origin of our brightest candidate. It completely sat- 
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Fig. 4. — Images of the z ~ 8 galaxy candidates identified in the CANDELS data. From left to right, these show: (1) a stack of all the 
optical ACS data (BViz), (2) WFC3/IR data in the filters Yios, (3) J125, and (4) ifi60- The images are 3 arcsec on a side, oriented North 
up, East left. All candidates are well detected in both J125 and H\qo, show a significant flux decrement to Y105, and show no flux in the 
optical. An unsmoothed version of these images are shown in the appendix in Figure [T2l 



isfies our selection criteria and our x%t limit such that 
it is most likely at z ~ 8, and so we just note the caveat 
raised by the clustering estimates. 

As can be seen from Figure HJ all z ~ 8 candidates 
are quite compact, but appear to be resolved. Their size 
distribution is consistent with the expectations based on 
extr apolation from lowe r redshift samples out to z ~ 7 
(e.g. Qes ch et al.l l2010ah . Only two sources show signs 
for elongated morphologies, indicative of disk structures. 
However, the surface brightness limit in the CANDELS 
data work against clearly detectin g such structures, c om- 
pared to e.g. the HUDF (see e.g. lOesch et al.ll2010aD . 

At the same time our paper was first submitted, the 
CANDELS team also published a c atalog of z ~ 8 ga laxy 
candidates from the DEEP area in lYan et al.l (pQlih . A 
significant fraction of their sources (45%) does not sat- 
isfy our strict optical non-detection criteria, with a few 
showing clear optical detections, thus disqualifying them 
for being at z > 7. All their candidates are discussed 
in detail in the appendix, where we also tabulate our 



measurements for these sources. 

3.5. Sample Contamination 

Thanks to our restrictive non-detection criteria we do 
not expect high contamination levels in our samples. 
Nevertheless, some residual contamination can not be ex- 
cluded. In particular, we briefly discuss several potential 
sources of contamination below: 

1. - Low Mass Stars: As shown in Figure [2j ultra- 
cool dwarf stars can exhibit similarly red colors as the 
potential high-redshift galaxies. Given the resolution of 
WFC3/IR, however, we can directly check whether any 
of our sources are unresolved. This test can reliably be 
done only for the brighter candidates, with well sampled 
profiles. Indeed, we did remove one source from our can- 
didate list with colors consistent with a T-dwarf which 
appeared to be point-like (at 03:32:25.33, -27:48:54.2). 
Furthermore, the contamination by dwarf stars is also 
expected to be low based on their obser ved surface densi- 
ties in high galactic latitude fields (e.g. Rv an et "aTl l2QQ5. 
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EoTTl). 

2. - Photometric Scatter: The most probable source 
of contamination in our sample comes from photometric 
scatter of sources with intrinsic colors similar to high- 
z galaxies. This has been quantified already in section 
13.31 To summarize, after applying the x%t < 3.0 cut, we 
only expect 1.0 ± 0.5 contaminants due to photometric 
scatter in the CANDELS Deep data set, and < 0.1 in the 
CANDELS Wide. Therefore, with our strict optical non- 
detection criteria, we limit the amount of contamination 
to - 10%. 

Similar conclusions are reached by investigating the 
xl pt distribution. As can be seen in the right panel of 

Figure [2j our sample includes five galaxies with x%t > 
0.9. Statistically, we expect only 22% of galaxies above 
such values, i.e. 3.5 sources. This would suggest again a 
contamination of ~ 10% in the full sample. 

3. - Spurious Sources: All our candidates are very 
well detected both in i^i6o,A# as well as in J125. The 
chance for spurious > 4.5a detections at the same loca- 
tion in both bands is negligibly small. Additionally, most 
of the sources are detected at lower significance in Y105, 
further reducing the chance of a spurious source. Spuri- 
ous sources are clearly not a concern for our sample. 

4. - Transients: Since the optical data has been 
taken a few years prior to the new WFC3/IR data, it is 
possible that supernovae that went off in the meantime 
are selected as infrared detections without any optical 
counterparts. However, since one of the specific goals of 
CANDELS-Deep is to search for such supernovae, the 
#160 an d J125 data acquisition has been distributed over 
several different epochs. So far, eight epochs have been 
acquired. This allows us to check directly whether any 
of our sources could potentially be a supernova. We 
therefore group the individual epochs in four bins sorted 
by exposure date and check whether all our sources are 
still visible in these sub-splits (in a #160 + J125 image). 
All CANDELS-Deep candidates are indeed detected in 
these images. For the two candidates in CANDELS- 
Wide data, this test is not as decisive, as the J125 and 
#160 data acquisition has only been split over two epochs. 
Furthermore, we stress again that these sources do not 
appear to be point-like. We therefore conclude that su- 
pernovae have not contaminated our sample. 



4. THE Z 



LF 



We will now use the z ~ 8 galaxy candidates identified 
in the previous section to derive constraints on the z ~ 8 
UV LF, and combine these new results with previous 
estimates to give the best available LF at z ~ 8. 

4.1. Selection Functions and Redshift Distributions 

To compute the LF, we first have to estimate the 
completeness, C( m), and redshift selectio n functions, 
SO, m). Following lOesch et al.l (|2007ll2009[ ). this is done 
by inserting artificial galaxies with varying magnitudes, 
profiles and sizes in the observational data and rerun- 
ning the source detection with the exact same setup as 
for the original catalogs. This is done for each of the 
fields individually. 

Our simulations are based on using real galaxies at 
lower redshift and scaling them to higher redshifts us- 
ing well-established evolutionary and cosmological re- 
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Fig. 5. — Top - The magnitude and redshift dependent selec- 
tion function of the CANDELS-Deep data as determined from our 
simulations. At bright magnitudes (#160, AB < 26 mag) > 50% of 
all simulated z = 7.3 — 8.9 galaxies are selected. At fainter mag- 
nitudes, this fraction is reduced due to photometric scatter in the 
color measurements and due to incompleteness. The selection func- 
tion for CANDELS- Wide is essentially identical, however, shifted 
by 0.7 mag due to the difference in the WFC3/IR data. Bottom - 
The redshift distribution of the CANDELS z ~ 8 galaxy candidates 
(blue thick line). The redshift distribution peaks at lower redshift 
than the selection function shown in the top panel due to the dim- 
ming of galaxies with redshift. The mean redshift of our sample is 
(z) = 7.9, with 80% of galaxies expe cted to lie at z = 7.2 — 8.7. 
The best-fit LF determined in Section ^31 was used for determining 
this redshift distribution. Also shown are the redshift distribution 
functions of the HUDF09 (gray solid) an d ERS (gray dashed) z ~ 8 
candidates from Bouwens et al. (2011b), which are very similar to 
the ones derived here for the CANDELS data. 



lationships, i.e . . usi ng the 'cloning' methodology of 
IBouwens et al.l ((2003). In particular, we 'clone' z ~ 4 
LBGs from the GOODS and HUDF fields to higher red- 
shifts. The images of these z ~ 4 sources are scaled to 
the desired input magnitude, and are stretched to ac- 
count for the difference in angular diameter distance, as 
well as a size scaling of (1 + z)~ l as observed for the Ly- 
man Break galaxy population across z ~ 3 — 7 (see e.g. 
jFerguson et all 120041 : iBouwens etHI 120041 : lOesch et al.l 
2010a). This procedure ensures that the distribution of 
morphologies and profiles of the simulated galaxies is 
as close to reality as possible and thus specifically ac- 
counts for SB-dimming effects, which result in loosing 
larger, resolved galaxies from the samples. The cloned 
galaxies are then inserted in the observed images with 
galaxy colors as expected for star-forming galaxies be- 
tween z = 6 and z = 9.5 (see also Figure [2]). The 
adopted colors are based on a UV continuum slope dis- 
tribution of j3 = — 2.5 ± 0.4 motivated by recent deter- 
minations of the UV continuum sl opes as a function of 
UV lu minosity at z > 6 (see e.g. [Bo uwens et al. 2009, 
2010b; Stanwav et al. 2005; Finkelstein et al. 2010, 2011; 
Wilkins et al. 2011b; Dunlop et al. 2012). As a cross- 
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Fig. 6. — The observed and the expected surface density of z ~ 8 
galaxies in GOODS-South. The two panels show the results for the 
CANDELS-Deep (upper) and CANDELS-Wide (lower) fields. The 
gray histograms show the observed surface density of candidates, 
while the two dark blue lines represent the expectations derived 
from our si mulations. The dashed l ine is derived from the best-fit 
z ~ 8 LF of Bouwens et al. (2011b). This somewhat overpredicts 
the observed number of sources at bright magnitudes found in our 
larger area survey. The solid line shows the expe cted surface den- 
sity using the new best-fit LF as derived in Section l4~4l The surface 
densities peak before the 5 sigma detection limits because of the 
reduction in the selection volume due to photometric scatter. 

check we have also tested that our simulation pipeline 
returns essentially equivalent results whe n adopting the- 
oretical Sersic galaxy pr ofiles (see also iBouwens et al.l 
I2Q11U [Qesch et alJl2Q12l ). 

From the simulation output, we compute the complete- 
ness as a function of observed #i60,A£ magnitude for 
each field, taking into account the scatter and bias be- 
tween input and output magnitudes. Additionally, we 
compute the selection probabilities as a function of red- 
shift and magnitude by measuring the fraction of sources 
that meet our selection criteria. The selection function 
and the corresponding redshift distribution function is 
shown in Figure [5] As can be seen, galaxies in our sam- 
ple are selected from z ~ 7.2 — 8.7, with a mean redshift 
of (z) = 7.9. 

We note that the expected redshift distribution has 
wings which extend to lower and higher redshift s, and 
we expect ~ 4 objects to lie at z < 7.5. The shape of 
the selection function has to be taken into account when 
modeling the observed number densities, particularly if 
there is strong evolution over the redshift range z = 7 — 
8. The redshift selection function is available from the 
authors upon request. 

4.2. Expected Surface Density of z ~ 8 Galaxies 

It is very instructive to compare the expected surface 
density of z ~ 8 galaxy candidates using the previous 
best-fit LF of IBouwens et all (|2011bf ) with the observed 
number of candidates in the CANDELS data. The ex- 
pected number of sources in a given magnitude bin mi 
can be estimated for any given LF, (j)(M), through: 

7Vf p = / dm f dz^S(m,z)C(m)<l)(M[m,z]) 

J Am J dz 

In Figure [6] we show the histograms of the observed 
surface density of sources in the different fields and com- 
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Fig. 7. The z ~ 8 UV LF from the total WFC3/IR data 
in GOODS-South. The stepwise LF from the CANDELS data 
alone is plotted as smal l bright bl u e circl es, while the previous 
determination from Bouwens et al. (2011b) using the ultra-deep 
HUDF09 and the ERS data is shown as light gray squares. The 
new CANDELS LF is clearly lower than the previous determination 
from the HUDF09+ERS at all magnitudes M uv < -20 by a factor 
~ 3 - 4. The best-fit LF from the HUDF09 and ERS data alone 
is shown as dashed gray line. The combined stepwise LF is shown 
as large blue circles, along with the best-fit as the solid blue line. 
The upper limit at the bright end corresponds to a la limit for a 
non-detection. This new LF was obtained by combining all the 75 
z ~ 8 candidates from the current CANDELS, HUDF09 and ERS 
data (see section [4~4|) . It is also listed in Table [3] For comparison, 
we also show the LF deter mination fro m the BORG survey as open, 
light blue squares Bradlev et al. (2012), which is in good agreement 
with our total LF determination. 

pare them to the expectation from different LFs. As can 
be seen, the expected source density peaks around 0.2 
arcmin -2 mag -1 in the CANDELS-Deep field. The peak 
occurs at #i60,ab = 27 mag, which is ~0.5 mag brighter 
than the formal 5<r limiting magnitude. This is mainly 
due to the reduction in the selection volume due to pho- 
tometric scatter, which causes us to lose sources from the 
color-color selection window near the limit. Note that we 
do also expect a small fraction of sources even below the 
formal magnitude limit due to scatter in the photomet- 
ric offsets between the aperture fluxes (from which the 
selection S/N is computed) and the total fluxes (used in 
Figure [6]) , as well as from variations in the depth of the 
#160 band data. 

From the integration of the surface density expected 
from the best-fit LF of IBouwens et all (|2Qllbf ). we find 
that we would expect to detect 22 z ~ 8 galaxies in 
the CANDELS-Deep field, and five in CANDELS-Wide. 
This is a factor 1.7 x higher than the 16 candidates we 
find, indicating that the bright end of the previous UV 
LF might have been estimated somewhat high due to 
the presence of a few very bright sources (#160, ab ~ 26 
mag) in a possible overdensity in the HU DFQ9-2 field (see 
also discussion in IBouwens et al.ll2011b[ ). The difference 
seems large but it is driven by just a few bright sources 
and so is consistent with small number statistics. 

4.3. New Constraints on the z ~ 8 LF 

We now use the new CANDELS z ~ 8 candidates 
to derive direct constraints on the UV LF, by comput- 
ing the step- wise LF in bins of absolute magnitudes. 
This is done using an approximation of the effective 
selection volume as a function of observed magnitude 
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Fig. 8. — The evolution of the Schechter function parameters and their uncertaintie s as a function of redshift. The contours show 68% 
and 95% of the likelihood at each redshift bin. The z < 8 contours are from[Bouwens et al. (2011b). Despite the large uncertainties, the 
z ~ 8 LF parameters are significantly different from z ~ 7. The difference is significant at > 99%. The characteristic magnitude is fainter 
by about 0.1 mag at z ~ 8, consistent with pure luminosity evolution, i.e. a dimming of M* to higher redshifts. The faint-end slope 
a is very steep, but it is still very uncertain at z ~ 8. Future deeper data will be required to improve on this (for reionization) crucial 
measurement of the slope a. 



V e fi(m) = / °° dz^S(z,m)C(rn). The LF is then given 

by <j>(Mi)dM = N? hB /Ves (ra»). This is shown in Fig- 
ure [TJ where we evaluated the LF in bins of 0.6 mag, as 
a compromise between luminosity resolution and S/N. 
The error bars include a 30% (40%) contribution from 
cosmic variance for CANDELS-Deep (CANDELS- Wide), 
w hich we estimate using t he cosmic variance c alculator 
of iTrenti fc Stiavellil (j20Q8h : see also iRobertsonl (|2010h . 

As expected from the low observed surface density of 
z ^ 8 candidates in the CANDELS data, the step-wise 
LF is significa ntly below the p reviou s best-fit Schechter 
function from IBouwens et all (j2011bh . At M uv = -21 
to —19 mag, we find values ~ 2 — 4x lower than from 
the HUDF09 and ERS data. As noted above, these are 
still consistent given the current and previous uncertain- 
ties. Nevertheless, the CANDELS data indicates that 
the characteristic magnitude might be fainter than pre- 
viously determined. We quantify this more precisely in 
the next section. 

Finally, we also compare our LF determination with 
the bright end co nstraints from the BORG survey 
(jBradlev et al.ll2Q12[ ) , which are based on 33 z ~ 8 galaxy 
candidates with J125 < 27.4 mag identified over an ef- 
fective search area of 274 arcmin 2 (i.e. about 1.7x the 
area of the HUDF09+ERS+C ANDELS) . As Figure 
shows, also the BORG z ~ 8 LF supports the indica- 
tion that the HUDF09 field is somewhat overdense in 
bright z ~ 8 galaxies. The LF from BORG is a factor 
~ 2 — 3 x lower than the step- wise LF determination from 
the HUDF09+ERS alone, and it is in excellent agreement 
with our new, total step- wise LF. 

4.4. Combination with Deeper Data: The z ~ 8 
Schechter Function 

Due to the small dynamic range in luminosities, the 
bright candidates identified in the CANDELS data are 
not sufficient to provide a good estimate of the overall 
shape of the UV LF alone. We therefore combine our 
new sources with all the z ~ 8 candidates from u ltra- 
deep field measurements from IBouwens et al.l ([201 lbl ) to 
update the Schechter function parameters of the z ~ 8 
UV LF. By doing so, we can now generate a LF using a 



TABLE 3 

Stepwise Determination of the z ~ 8 UV LF Based 
on CANDELS, HUDF09, and ERS Candidates 



M uv [mag] 


</>* [10- 3 Mpc- 3 mag _1 ] 


-21.94 


< 0.008* 


-21.34 


0.005 ± 0.006 


-20.74 


0.046 ± 0.020 


-20.14 


0.130 ± 0.045 


-19.54 


0.339 ± 0.118 


-18.94 


1.03 ± 0.35 


-18.34 


1.56 ± 0.72 


-17.74 


4.52 ± 2.07 



* la upper limit for a non-detection. 



total of 75 z ~ 8 candidates identified over 148 arcmin 2 , 
spanning i^i60,AB ~ 25.7 — 29.5 mag. The stepwise de- 
termination of this combined sample is listed in Table [3] 
and shown in Figure [71 

The Schechter function parameters are derived by max- 
imizing the Poissonian likelihood for observing 7V obs 
sources in a given magnitude bin when 7V exp are expected 
to be seen based on a given UV LF. We thus maximize C 

= Ylj Ui P ( N j h i S ' N j7)' where 3 runs over a11 fields ' and 
i runs over the different magnitude bins, and P is the 
Poissonian probability. The expected number of sou rces 
are computed according to the equation from section H~2l 
For the HUDF09 and ERS fields we adopt t he selection 
volumes estimated in IBouwens et al.l (|2011bf l . The com- 
bination of these older results with our new sample is 
appropriate since both are based on essentially identical 
simulations and inputs, i.e., our estimat es of the CAN- 
DELS s election functions match those of IBouwens et al.l 
(f20lTbh . 

The best-fit parameters are determined by a grid 
search over Schechter function parameters, maximizing 
the combined likelihood C. The best-fit solutions are: 
log(0* [Mpc^mag" 1 ]) = -3.30™, M* = -20.04^;^ 

mag, and a = — 2.06^q'28« These parameters are consis- 
tent with the previous de termination from the HUDF09 
and ERS data alone (jBouwens et al.ll2011b[ ). and, despite 
a somewhat larger sample size, the uncertainties are not 
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TABLE 4 

Comparison of z ~ 8 LF Determinations in the Literature 

Reference log</>* [Mpc _3 mag _1 ] ^\jv [ ma §] a 

This Work -3.30to 46 -20.04^;^ -2.06^0 28 

Bradley et al. (2012) - 3 - 37 ±o'.2i - 20 - 26 ±o'.34 - L98 ±o.22 

Bouwens et al. (^Ollb) -3.23±Q 27 -20.10 ±0.52 -1.91 ± 0.32 

Lorenzoni et al. (2011) -3.0 ' -19.5 -1.7 (fixed) 

Trenti et al"7201l£) -3.4 (fixed) -20.2 ± 0.3 -2.0 (fixed) 

McLure et al. (2010) -3.46 -20.04 (fixed) -1.71 (fixed) 

Bouwens et al.72010a) -2.96 (fixed) -19.5 ± 0.3 -1.74 (fixed) 



markedly reduced. However, the best-fit characteristic 
magnitude is fainter by ~ 0.1 mag, which is due to the 
lower number of detected sources in the CANDELS data 
than the number that was expected from the previous 
LF determination. 

Note also that the fai nt-end slope is steepe r than 
a = —2, as discussed by [Bouwens et all (|2Q12a[ ). This 
leads to a formally divergent luminosity density. How- 
ever, galaxies are not expected to be formed below a 
given luminosity due to inefficient cooling in low mass 
halos and feedback effects (at masses that correspond to 
Muv about —10 to —11), and so the luminosity density 
converges. Nevertheless, such steep slopes have impor- 
tant consequences f or reionization by gala xies as pointed 
out previously (e.g. IBouwens et al. 2012a). 

At the bright end (M uv < -20), our best-fit LF is a 
factor ~ 1.7 x lower than the previous determination of 
IBouwens etHI (|20Ub|), who used only the HU DF09 and 
ERS data. However, IBouwens et all (j2011bh only find 
such a higher surface density for z ~ 8 sources in the 
HUDF09, while the number counts over the ERS are ac- 
tually in excellent agreement with what we find over the 
CANDELS field. Therefore, all the wide-area data over 
the CDF-South GOODS yield approximately the same 
surface density of bright z ~ 8 galaxies. Nevertheless, 
we stress again that both the step-wise and the best-fit 
LFs are consis t ent wit h the previous determination of 
IBouwens et al.l (j2011bl ). given the still small number of 
sources at the bright end. 

In Table [4j we compare our new best-fit LF parame- 
ters w ith other, previous de t erminations from the liter- 
ature (IBouwens et al l 2010a, 2011b; M cLure et aHl2010l: 
Lorenzoni et al.ll20lH iTrenti et alJl2011bt iBradlev et all 
2012f ). Within the current measurement uncertainties 
these are all consistent with each other. 

Overall, it is very reassuring that different groups ar- 
rive at similar results given the variety of approaches. 
Nonetheless, it is worthwhi l e rem arking t hat with the 
except ion of IBouwens et al.l (j2Qllbh and IBradlev et al.l 
(2012) the previous determinations were not based on 
large enough data sets such that all three Schechter func- 
tion parameters could reliably be fit simultaneously as 
done here. 

4.5. Evolution of the LBG Population at 4 < z < 8 

A key diagnostic in studying the build-up of galaxies is 
how the UV LF evolves with redshift, as this is directly 
related to the distribution of SFRs in galaxies. In Figure 
El we show the error contours of the Schechter function 
parameters from current HST data of LBGs at different 
redshift s. The error contours for the 4 < z < 7 samples 



are taken from IBouwens et al.l (|2007l 1201 lbh . while the 
z ~ 8 contours correspond to the likelihood contours of 
our Schechter function fit at z ~ 8 from all available data, 
i.e. including the CANDELS GOODS-South candidates 
(Section @~4]). 

The Schechter function parameters show significant 
evolution from z ~ 7 to z ~ 8. In particular, the combi- 
nation of the characteristic magnitude and number den- 
sity is evolving at > 99% significance. On the other hand, 
the combined constraint on the faint-end slope and the 
characteristic luminosity is still quite weak (right panel 
of Figure [8]) . Unfortunately, the faint-end slope is largely 
unconstrained. This is very unfortunate since it is one 
of the most important parameters for assessing the con- 
tribution of galaxies to reionization. The flux density 
of ionizing photons from galaxies is extremely sensitive 
to the faint-end slope a , especially when the slope is as 
steep as a ~ — 2 (e.g. IBouwens et al.l l2012al ). A con- 
tinued effort in the future will thus be to constrain this 
parameter better with deeper WFC3/IR imaging. 

Overall, the UV LF constraints of LBGs are con- 
sistent with pure luminosity evolution from z ~ 8 
to z ~ 4. This can easily be a chieved by grow- 
ing SFRs of in dividu al galaxies (e.g. Stark et al.l 120091: 
iFinlator et all l2QllU _ jPapovich et all 120111 : iSmit et all 
120121 : Uaacks et all 1201 2( 1. The best-fit evolution of M* 
of LBGs at z > 3 follows: 

M*(z) = -20.98(±0.04) + 0.31(±0.03) x (z - 3.8). 

This is shown in Figure [9j where we plot the evolution of 
the characteristic cut-off luminosity of the UV LF with 
redshift. After correcting for dust extinction, the SFR of 
an L* galaxy thus grows by almost an order of magnitude 
from z ~ 8 to z ~ 4. 

The value of M* as measured at z ~ 8 is essentially 
equal to the cut-off UV luminosity at z ~ 1.5-2, just 
shortly after the peak of the cosmic SFR density. How- 
ever, the SFR of an L* galaxy is nevertheless larger at 
z ~ 1 .5-2 than at z ~ 8 due to the larger dust obscuration 
^ gJReddy et al.ll2010l : IBouwens et aLll2Q12bt ISmit et al.l 

120m 

In Figure [10] we show how the (dust-corrected) SFR 
density evolves across redshift when including our new 
UV LF parameters and the updated dust corr ections 
for LBGs at z > 4 from IBouwens et all (|2Q12bh . The 
SFR density includes all galaxies down to a fixed flux 
limit of Muv = —17.7, which is the current detection 
limit at z ~ 8. Our updated SFR density at z ~ 8 
is only marginally lower than the previous measure- 
ment of IBouwens et all (|2011bf ) by 0.07 dex, i.e. we de- 
rive log psfr = -2.32 ± 0.12 Moyr^Mpc- 3 . This is 
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Fig. 9. — Evolution of the characteristic luminosity of the UV LF 
across redshift. The measurements at z > 3.5 shown as dark red 
circles are based on LBG selections in the deepest HST images, 
including the HUDF and H UDF09 fields. These are taken from 
IBouwens et al.l (2007, 2011b). In general, determinations by other 
authors are in good agreement with these measurements. The red 
square is the best-fit value for M* at z ~ 8 when combining the 
HUDF09+ERS data w ith o ur new measurements from the CAN- 
DELS data (see section |4.4|) . The lower redshift measurements are 
a sele ction of UV LF parameters determined by Reddv & Steidel 
(2009, gray diamon ds) at z ~ 2 - -3, and Oes ch et al.l (1201 Obi . 
black open squares), Arnouts et al. (2005|, black open circles), and 
ICucciati et al.l (2012, black crosses) at z < 2. The dashed black 
line and gray shaded area correspond to the best-fit evolution of 
M* as a function of redshift at z > 4. 
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Fig. 10. — The evolution of the star-formation rate density with 
redshift. The measurements are based on the integration of the UV 
LFs to Mjjy = —17 .7 and converting to SFRs using the relation of 
Madau et al. (1998). The upper, red shaded area is obtained after 
correcting the observed UV luminosity densities for dust extinction 
based o n the UV continuum slope distributions of Bouwens et al. 
( 2012b|). The plot symbols are the same as in Figure [9] The SFR 
density estimate at z ~ 10 is based on the d etection of one galaxy 
candidate in the HUDF (jOesch et al.l [2011 ). The rapid (but still 
uncertain) increase in the cosmic SFR density from z ~ 10 to z ~ 8 
will be refined with upcoming WFC3/IR imaging over the HUDF. 

based on the luminosity density of log pL = 25. 58 ±0.12 
erg s~ 1 Hz~ 1 Mpc ~ 3 and the conversio n of the UV lumi- 
nosity to SFR by lMadau et al.l (|1998[ ). 

As can be appreciated from Figure [TOj the SFR density 
grows rather dramatic ally by more than an order of mag- 
nitude from z ~ 10 ( Oe sch et al.l 120121 : IBouwens et al.l 
l2011af ) to our new determination at z ~ 8. After that it 
evolves very steadily, growing by another ~ 1.5 dex from 
z ~ 8 to its peak at z ~ 2.5. The rapid growth at z > 8 
is very intriguing but still very uncertain. The change 
at z > 8 may well be refined with planned deep F140W 
data over the HUDF. 



5. SUMMARY 

In this paper, we analyzed newly acquired WFC3/IR 
data over the GOODS-South field as part of the CAN- 
DELS MCT program to derive new constraints on the 
UV LF at z ~ 8, about 600-750 Myr after the Big Bang. 
Unlike at z ~ 7, where the brightest galaxies can be 
detected from the ground, z ~ 8 galaxies are too faint 
to be detected with any reliability with current ground- 
based instrumentation. Therefore, the bright end of the 
z ~ 8 UV LF has remained relatively uncertain since it 
has been based on WFC3/IR data that was only avail- 
able over a relatively small area, the very deep HUDF09 
fields. 

Galaxies at z ~ 8 are identified using the Lyman Break 
technique on optical and near-IR data, and particularly 
the Y105 band, leading to them being identified as 'Y105- 
dropouts'. The CANDELS fields over GOODS are ex- 
tremely valuable data sets, as they combine imaging in 
three WFC3/IR filters, and in addition also have rela- 
tively deep ancillary optical data from HST ACS, which 
allows for such z ~ 8 galaxy selections. Our selection cri- 
teria are outlined in section [3] and were chosen to match 
prior measurements in the ultra-deep HUDF09 dataset. 
They select galaxies at z ~ 7.2 — 8.7 with a mean redshift 
(z) = 7.9. 

In the full search area of ~ 95 arcmin 2 , we identified 16 
new z ~ 8 galaxy candidates with i^i60,A£ magnitudes in 
the range 25.7 — 27.9 mag. These are presented in Figure 
H] and Table El The z ~ 8 candidates were selected using 
very strict optical non-detection requirements, incl udin g 
a measurement of the optical xl p t ^ ux ( see Sectiorj3J3]) • 
This allows us to reduce the expected contamination due 
to photometric scatter and due to the limited depth of 
the optical data relative to the WFC3/IR imaging by a 
factor ~ 3 — 4. Our sample is thus expected to show 
relatively low levels of contamination, which we estimate 
to be - 10%. 

Interestingly, the observed surface density of z ~ 8 
galaxies is ~ 1.7x lowe r than expected from t he previ- 
ous best-fit z ~ 8 LF of IBouwens et al.l ([201 lb[ ). though 
the difference is consistent with the statistics from the 
small numbers of objects at bright magnitudes in the 
previous search. This previous estimate was based on 
the much smaller area data from the HUDF09 fields and 
the ERS data. By using the current CANDELS imag- 
ing we triple the search volume for luminous galaxies at 
Muv < —19.5 mag in the CDFS relative to our previous 
analysis, thus reducing the potential biases introduced 
by cosmic variance. 

We combine our new CANDELS z ~ 8 candidates with 
our previous candidates from the much deeper HUDF09 
and ERS fields to derive the best possible measurement 
of the UV LF. The best-fit z ~ 8 UV LF we derive in 
this way is con s istent w ith the previous estimates from 
IBouwens et al.l (|2011bl ). However, the best-fit charac- 
teristic magnitude is fainter by ~ 0.1 mag (M* =s = 
—20.04 ±0.46 mag). Despite the larger area probed here, 
the uncertainties on the Schechter function parameters 
are still significant. Nonetheless, these new results are 
contributing to our growing understanding of the evolu- 
tion of key parameters in the luminosity function from 
z ~ 8 to z ~ 2, and of the star formation rate density 
from z ~ 10 to z ~ 2. 
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A combination of future wide area data, as well as 
deeper WFC3/IR imaging will be necessary to further 
improve on these very important measurements that are 
crucial for refining our estimates of the role of UV pho- 
tons from galaxies in the reionization of the universe be- 
fore the advent of JWST. 

We are grateful to Massimo Stiavelli, Kristian Finlator, 
Steve Finkelstein and Naveen Reddy for helpful discus- 
sions related to this work. We also thank the CANDELS 



collaboration for their big efforts in planning and ob- 
taining such an exquisite data set available to the com- 
munity. Support for this work was provided by NASA 
through Hubble Fellowship grant HF-51278.01. This 
work has further been supported by NASA grant HST- 
GO-1 1563.01. This research has benefited from the SpeX 
Prism Spectral Libraries, maintained by Adam Burgasser 
at |http: / /www .browndwarfs.org/spexpr ism[ 
Facilities: HST( ACS / WFC3) , Spitzer(lRAC). 
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APPENDIX 

A. COMPARISON TO YAN ET AL. (2011B) CATALOG 

At the same t ime our paper was submitted, the CANDELS team submitted a very similar z ~ 8 galaxy search in 
lYan et all (|2QlTf ). That paper is still under review and is likely to change further (Yan 2012, private communication). 
Below we explicitly compare onl y to version 2 of their manuscript, which is publicly available on arXiv at the time 
of publication of our paper. The lYan et all (j2011f ) analysis is not based on the complete GOODS-South CANDELS 
data set (they include 83% of the final data in J125 and Hiqq that was used in our analysis). Additionally, they only 
analyze the CANDELS-DEEP area (see Fig [I]), and they only include 8-epoch mosaics of the Jgi4 image. 

For comparison, our analysis is based on all the isi4 data available over the GOODS-South, which was obtained from 
a large number of additional programs (ERS, HUDF09, CANDELS SNe follow-up, 3D-HST, and the UVUDF). As 
noted in the main text, these data reach almost 1 mag deeper than the GOODS 2775 data, and are thus an extremely 
valuable addition to check for low-redshift contaminants. Additionally, we include a deeper reduction of the other ACS 
images including additional data taken after the original GOODS-South program, which reaches to ~ 0.1 — 0.2 mag 
deeper compared to the GOO DS v2.0 dat a. 

The revised v2 manuscript of I Yan et al.l (|2011[ ) lists a total of 16 different Yios-dropout candidates, which are split in 
two different samples, one using AUTO fluxes and one usi ng ISO fluxes for computing colors. Out of these 16 sources, 
only 5 are in common with our sample. lYan et al] (|2011f ) use somewhat more restrictive color criteria than what we 
adopt here (Y105 — J125 > 0.8, instead of > 0.45 as used in our analysis). Therefore, we would have expected to select 
all their sources in our catalog. 

Upon inspection of their sources, we found the main reason that I Yan et al.l (j2011f ) z ~ 8 candidates did not appear in 
our catalog is due to a large fraction of their candidates showing non-negligible flux in the optical ACS data such that 
they do not satisfy our strict optical non-detection criterion x%t < 3. This applied to 7 out of their 16 candidates. 
lYan et aTl (j2011l ) state that they perform a two-stage source selection, which includes a visual inspection of their 
candidates using all the data. However, some sources show detections in our reduction of the optical ACS data, 
disqualifying these sources being at z > 7. Furthermore, two sources are detected in the HUDF ACS data. Since these 
data were used in their visual inspection, we do not have an explanation for this discrepancy. 

In addition to these 7 sources with optical flux, their list also includes two sources which are not significantly detected 
(< 5<j) in our i^i6o-band image, as well as two sources for which we measure Y — J colors that are too blue, although 
they are consistent with being y-dropouts at the la level. 

A detailed summary of all our measurements and analysis of the Yan et al.l (|2011f ) sample is provided in Table [5j 
and an example of optical stamps for two sources are show n in Figure TT1 

Given this disagreement, it is encouraging, however, that I Yan et al.l (j2011f ) do, in fact, include 5 out of the 6 sources 
with S/N(Hiqo) > 7 from our sample which satisfy their stricter color criterion of Y — J > 0.8 (see Table [2]). This 
indicates that our selection is quite robust, at least for higher-significance sources, if optically detected interlopers are 
excluded properly, as done for our catalog. 

We stress again that making full use of all the information in the optical data is extremely important for a reliable 
LBG selection. We did our best possible effort to do this by including a limit in the optical Xo P t measurement and by 
analyzing all the available optical ACS data that were taken over this field. Furthermore, we note that as long as the 
contamination fractions and detection efficiencies are modeled self-consistently, using shallower data should not result 
in a different estimate of the final LF (within the errors). 
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TABLE 5 

Comparison to Source List of Yan et al. (ApJ submitted, arXiv:1112.6406v2) 



Von TTt 




a 

















J\2b 






Q /"W -r-r 


«^125 — ^160 


*105 — ^125 


Note 


ISO.085 


03 


32 


49 


94 


-27 


48 


18 


1 


25 


8 


1 

± 


0.1 


(25 


6) 


22 


0.02 ± 0.08 (-0.1) 


1.00±0.11 (1.20) 


OK, in our catalog 


ISCL164 


03 


32 


41 


42 


-27 


44 


37 


8 


26 


1 


± 


0.1 


(26 


0) 


17 


-0.00 ± 0.10 (0.0) 


0.82±0.14 (0.90) 


y2 = 3 3 > 3 


ISO_157 


03 


32 


42 


88 


-27 


45 


04 


3 


26 


g 


± 


0.1 


(26 


5) 


10 


U.U 1 _1_ u. 10 IU.1 1 


n 7«4-0 27 CO QO") 


v 2 — 8 7 > 3 
X op t — °- ' ^ ° 

OK, in our catalog 


TQM 071 


Uo 


32 


20 




-A I 




00 





Z t 


1 
1 


± 


0.2 


( 9£i 

\Zv 


y ) 


o.z 


U.U/ zt U.zo (,U.1J 


1 QO-Lfl /IS 
l.oZztU.4o (^>1.<UJ 


ISO.082 


03 


32 


14 


13 


-27 


48 


28 


9 


26 


9 


± 


0.1 


(27 


2) 


7.1 


-0.47 ± 0.24 (0.0) 


>2.1 (>1.10) 


S/N(W) = 2.1 


ISO_078 


03 


32 


41 


65 


-27 


48 


34 


5 


27 


4 


± 


0.2 


(27 


2) 


5.0 


-0.01 ± 0.26 (0.2) 


>2.3 (>1.70) 


X2pt = 4.6 > 3 a 


ISO_011 


03 


32 


14 


47 


-27 


51 


48 


5 


27 


4 


± 


0.3 


(27 


2) 


3.6 


-0.35 ± 0.39 (0.0) 


>1.7 (1.20) 


too low S/N in ifi6o 


ISO_158 


03 


32 


47 


95 


-27 


44 


50 


4 


27 


5 


± 


0.2 


(27 


3) 


7.6 


0.34 ± 0.24 (0.2) 


0.39±0.30 (0.80) 


too blue Y — J 


ISO_063 


03 


32 


35 


00 


-27 


49 


21 


6 


27 


3 


± 


0.2 


(27 


6) 


7.1 


0.10 ± 0.22 (0.1) 


>2.3 (1.50) 


OK, in our catalog 


ISO.017 


03 


32 


18 


09 


-27 


51 


18 


5 


27 


9 


± 


0.2 


(27 


6) 


4.0 


-0.04 ± 0.31 (-0.4) 


>1.5 (1.30) 


X%t = 3-01 > 3 


ISO.160 


03 


32 


46 


11 


-27 


44 


48 





27 


8 


± 


0.4 


(27 


9) 


3.8 


-0.07 ± 0.33 (0.1) 


0.20±0.30 (0.90) 


S/N z = 2.8, S/N; = 3.1 b 


ISO.008 


03 


32 


16 


91 


-27 


52 


01 


9 


27 


7 


± 


0.3 


(28 


0) 


2.7 


-0.75 ± 0.57 (-0.3) 


>1.5 (0.90) 


too low S/N in ifi60 


AUTO_212 


03 


32 


20 


96 


-27 


51 


37 


1 


26 


4 


± 


0.1 


(26 


4) 


7.8 


-0.19 ± 0.16 (0.1) 


1.15±0.26 (1.00) 


OK, in our catalog 


AUTO_368 


03 


32 


34 


50 


-27 


46 


03 


5 


27 


1 


± 


0.2 


(27 


2) 


7.5 


0.30 ± 0.19 (0.2) 


0.61±0.33 (0.80) 


S/N z = 2.6, S/N; = 2.1 c 


AUTO_204 


03 


32 


18 


18 


-27 


52 


45 


6 


27 


3 


± 


0.2 


(27 


4 ) 


7.0 


-0.16 ± 0.23 (0.2) 


1.06±0.44 (0.90) 


OK, in our catalog 


AUTO.094 


03 


32 


40 


67 


-27 


45 


11 


6 


27 


5 


± 


0.2 


(27 


4) 


4.1 


-0.24 ± 0.34 (0.0) 


0.28±0.30 (0.90) 


too blue Y — J 



* 

The values in this table are our own measurements. For comparison, we show the measurements from Yan et al. in parentheses. Limits are la. 
a Source ISO.078 is detected in several optical filters of the HUDF ACS data. 

k Source ISO_160 is detected in all optical bands of our reduction of the GOODS data at more than 1.5<r. 
c Source AUTO_368 is additionally detected in the V 606 data of the HUDF. 
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Fig. 11. — Two examples of clear optical detections in the Yan et al. (2011) sample (ISO_160, and ISO_164). Both sources show a 
detection in Isi4, and show hints of flux in other bands as well. The increased depth of the /si4 data is clearly an extremely valuable 
addition for removing contaminants. Note that both these sources were flagged as contaminants in our catalog based on our Xopt < 3 
criterion, showing the power of using this measurement. 



16 



Oesch et al. 
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Fig. 12. — Same as Figure 0] but without the slight smoothing of the stamps. Smoothing more accurately displays a larger dynamic 
range in the images and thus shows the reliability of the sources more fairly. 



